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1.0  INTRODUCTION 


At  the  present  time,  there  is  no  generally  acceptable  model  of  sea  ice 
based  upon  the  physics  of  the  material.  Two  models  have  been  proposed  in  the 
past,  and  both  models  fall  short  of  adequately  explaining  the  strength  or 
elasticity  of  sea  ice. 


2.0  MODELS  OF  SEA  ICE 


2.1  STRENGTH  MODEL 

The  first  model  was  developed  by  Anderson  and  Weeks  (1).  This  model  was 
based  upon  the  experimental  observations  of  Weeks  and  Anderson  {2}.  From  their 
observations  they  deduced  that  the  salinity  and  temperature  were  the  controlling 
parameters  for  a  sea  ice  strength  model.  Their  experimental  observations  showed 
them  that  brine  was  trapped  between  crystal  boundaries  and  subcrystalline  plate¬ 
lets  in  vertical  layers.  With  the  nucleation  and  formation  of  sea  ice,  and  then 
with  further  decrease  in  temperature,  the  brine  layers  start  to  decrease  in 
cross-sectional  area.  Initially,  with  decrease  in  temperature,  the  brine  layers 
start  to  neck,  then  form  vertical  cylinders,  and  finally  separate  into  elliptical 
cylinders  and  circular  cylinders.  This  sequence  of  events  with  the  sea  ice  brine 
lavers  is  shown  in  Figure  1. 


Anderson  and  Weeks  reasoned  that  if 


-  d/1 

o 


(1) 


where  l  is  brine  in  a  unit  area,  d  is  the  diameter  of  the  brine  cvlinders,  and 
e 

1  is  the  separation  of  cylinders  in  the  failure  plane  (see  Figure  1) ,  then  the 
strength  of  sea  ice  a  will  be  a  function  having  the  following  form 

3  ^  o  (1  -  8  )  (2) 

P  e 


where  a 

P 


is  the  strength  of  pure  ice. 
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Figure  1.  Schematic  diagrams  showing  how  a  brine  layer  splits  into 
cylinders.  (From:  Anderson  and  Weeks',  see  reference  ;1}) 


The  stress  concentration  factor  k  is  defined  as  the  maximum  stress  to  the 
average  stress  in  the  minimum  cross  section  (3,  4}. 


k  *  .  P 
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(3) 


or  for  the  sea  ice  strength  we  will  have 

J  ■  o  (1  -  8  )/k  .  (4) 

P  e 

Next  Anderson  and  Weeks  considered  round  cylinders  and  the  geometrical  properties 
of  sea  ice  to  develop  the  following  relation  for  the  brine  volume  S,  i.e.,  the 
volume  of  the  brine  to  the  total  volume 

2  ■  ir8^/4a  (5) 

e 
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where  W  ■  1  a  and  W  is  the  platelet  width.  Hence,  from  Eqs.  (3)  and  (4)  the 
relation  for  sea  ice  strength  can  be  written  as 


(6) 


where  a  is  a  function  of  brine  pocket  separation.  This  approach  considered  the 
sea  ice  strength  to  be  a  function  only  of  brine  volume,  brine  pocket  separation, 
pure  water  ice  strength,  and  a  nondefined  stress  concentration  factor.  Brown  ' o) 
shows  a  comparison  of  his  experimental  data,  shown  in  Figure  2,  with  the  Anderson 
and  Weeks'  model.  In  Figure  2,  a is  the  flexural  strength  of  pure  ice  and  k  is 
the  stress  concentration  factor.  The  comparison  shows  there  is  considerable  room 
for  improvement. 


2.2  ELASTICITY  AND  STRENGTH  MODEL 

Following  Anderson  and  Weeks'  work.  Brown  {5}  considered  a  perforated  plate 
model  as  is  shown  in  Figure  3  for  a  sea  ice  model.  The  plate  is  infinite  in 
extent,  and  the  extensions  in  x  and  y  are  uniform,  then  any  stresses  produced 
around  any  one  hole  are  harmonic  with  respect  to  those  of  the  other  holes.  For 


FLEXURAL  STRENGTH  (KG/CM2) 


Figure  2.  Comparison  of  experimental  flexural  strength 
versus  brine  content  in  low-temperature  layer  with  the 
Anderson  and  Weeks'  theoretical  model. 
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Figure  3.  Plan  view  of  a  section  of  a  perforated  plate  with 
a  square  array  of  holes.  [From:  Bailey,  R.  and 
R.  Hicks,  J.  Mech.  Eng.  Sci.,  2(2),  143  (I960).] 

analysis,  one  need  consider  only  one  square  element  with  the  coordinates  origi¬ 
nating  at  the  center  of  the  hole.  In  approximating  a  perforated  plate,  it  is 
assumed  the  internal  properties  of  sea  ice  have  brine  cylinders  in  which  the 
diameters  are  very  small  as  compared  to  the  lengths  of  the  cylinders.  In  addi¬ 
tion,  the  lengths  of  the  cylinders  are  normal  to  the  plane  of  the  sea  ice  sheet 
and  the  cylinders  at  different  levels  of  the  ice  sheet  add  to  become  perfora¬ 
tions  of  uniform  extent  throughout  the  ice  sheet.  Hence,  we  assume  Bailey  and 
Hicks  {6}  perforated  plate  theory  is  valid  for  sea  ice.  From  Bailey  and  Hicks' 


theory,  let’s  consider  the  perforated  plate  as  shown  in  Figure  3  to  be  in 
tension.  Then  from  Bailey  and  Hicks'  theory,  we  have  the  following  for  Young's 
modulus  £  and  Poisson's  ratio  v: 


2r  Z 
12 


:1  +  Z2 


r  -  r 

2i _ 21 

Zl  +  I2 


where  E'  refers  to  Young's  modulus  of  a  material  without  any  holes  (pure  ice  in 
this  case) ,  and  where 
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where  fell  symbols  are  defined  in  the  appendix) 
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and 
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where  the  primes  refer  to  a  material  without  holes.  Assume,  as  Anderson  and 
Weeks  did,  that  sea  ice  failure  in  tension  occurs  across  a  row  of  brine  cylinders 
that  are  in  the  minimum  cross-sectional  plane  (the  ligament  between  cylinders) 
passing  through  the  circular  cylinders  normal  to  the  applied  stress.  If  the 
failure  stress  is  the  average  stress  across  the  ligaments  at  failure  and  if  that 
failure  occurs  with  an  elastic  break  (i.e.,  on  the  linear  portion  of  the  defor¬ 
mation  curve),  then  from  Bailey  and  Hicks  (6)  we  have  the  ratio  of  the  average 
stress  across  the  perforated  plate  ~  to  the  stress  of  a  plate  without  any  holes 
a'  as  given  by 


rr  =  k 


i/9VEVc 

2\  Z1 


9e./EU  / c ’ 
2  x 


in  which  k  is  the  average  stress  concentration  across  the  ligament  and,  with  an 
angle  of  9  =  tt/2,  90  becomes 


+  (n  +  2) (n  +  1)B 


+  (n  -  2)  (n  -  1)  B 


r.(f)  +  n(n  +  1)An(f) 

/a\n| 

n(r J  “s  n9 


and  wherein  99^/(EUx/c)  has  values  of  n  =  0,  4,  8,  •  *  •  and  902/(ElWc)  has 
values  of  n  *  2,  6,  10,  •  •  •. 


Bailey  and  Hicks  have  calculated  the  coefficients  for  Eqs.  11,  12,  and  17 
and  have  graphed  most  of  the  functions.  The  bars  over  several  of  the  functions 
indicate  averaging  which  has  been  done  along  the  boundaries  and  ligaments  of  the 
unit  plate. 


From  Eqs.  16  and  17  it  is  possible  to  determine  the  relative  plate  wave 
velocities  and  the  relative  flexural  strengths  with  the  conditions  assumed  here. 
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A  comparison  of  the  perforated  plate  theory  with  experimental  field  data  of 
relative  plate  wave  velocity  (Vp/V^) ,  for  elasticity,  is  shown  in  Figure  4.  it 
should  be  noted  that  the  ratio  c/a  is  the  ratio  of  one-half  the  unit  cell  side  to 
the  radius  of  the  brine  cell.  Figure  4  shows  the  comparison  between  theory  and 
experimental  data.  While  the  comparison  is  not  good,  it  does  show  that  elasticity 
is  a  function  of  brine  content. 

A  comparison  of  the  perforated  plate  theory  with  the  flexural  strength  of 
sea  ice  is  shown  in  Figure  5.  Again,  the  comparison  indicates  that  flexural 
strength  is  a  function  of  brine  content;  however,  there  must  be  other  factors 
influencing  the  flexural  strength. 

A  comparison  of  theory  and  experimental  values  of  plate  wave  velocitv  with 
flexural  strength  is  shown  in  Figure  6.  The  comparison  between  theory  and  exper¬ 
iment  appears  to  be  relatively  good  for  geophysical  data. 


RELATIVE  PLATE  VELOCITY  (Vp/V'p) 

Figure  4.  Comparison  of  experimental  data  with  perforated  plate 
theory  for  brine  content  versus  the  relative  plate 
wave  velocity,  with  V'  *  3,410  m/sec. 


e/a =1.2 


A  -Wales,  Alaaka:  1958 
•  -Walaa,  Alaska:  1959 
■  -Thule,  Or  sent  and:  1961 


RELATIVE  FLEXURAL  STRENGTH  (cr/cr') 


Figure  5.  Comparison  of  experimental  data  with  perforated  plate 
theory  for  brine  content  versus  the  inverse  of  the  relative 
flexural  strength,  with  o?  *  8.0k  kg/cm2. 


Figure  6.  A  comparison  of  plate  wave  velocity  versus  flexural  strength 
from  theory  of  perforated  plates,  with  the  experimental 
values,  letting  V'  ■  3,410  m/sec  and  o’  *  8.0k  kg/cm^. 


3.0  SUMMARY 


This  review  of  theoretical  models  of  elasticity  and  flexural  strength  of 
sea  ice  is  to  establish  a  baseline  for  future  efforts  in  experimental  and  theo¬ 
retical  work  in  these  areas.  For  example,  brine  content  is  not  the  total  porosity 
of  sea  ice.  Total  porosity  of  sea  ice  is  brine  plus  gas  content.  Other  factors 
which  might  influence  the  models  and  should  be  considered  in  any  theoretical 
models  for  sea  ice  strength  and  elasticity  are: 

1.  Visco-elastic  and  plastic  effects  along  the  failure  plane  as  well  as 
slippage. 

2.  The  distribution  and  geometry  of  the  pores. 

3.  Gas  porosity  must  be  taken  into  consideration. 

4.  Consider  other  configurations  for  theoretical  models. 

5.  Surface  tension  effects  on  the  cohesion  of  the  ice  platelets. 

6.  Effects  of  solid  salts  on  cold  ice. 
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APPENDIX 

DEFINITION  OF  SYMBOLS 


A  ,  B  ,  etc 
n  n 

a 

c 

E 

EE 
r,  e 

r\  r\  r\ 

rr,  99,  r0 


x,  y 

a  n 
XX,  vy ,  xy 

v 

v  v 


Constants  in  stress  function 
Radius  of  hole 
Pitch  of  hole 
Young's  modulus 
Effective  Young's  modulus 
Polar  coordinates 

Stress  components  in  polar  coordinates 

Displacements  in  cartesian  coordinates  for  a 
plate  with  a  diagonal  system  of  holes 

Cartesian  coordinates 

Stress  components  in  cartesian  coordinates 
Poisson's  ratio 

Numerical  integrations  as  defined 
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